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ABSTRACT 

We report on new results from simultaneous, dual frequency, single pulse 
observations of PSR B0329+54 using the Giant Metrewave Radio Telescope. We 
find that the longitude separation of subpulses at two different frequencies (238 
and 612 MHz) is less than that for the corresponding components in the average 
profile. A similar behaviour has been noticed before in a number of pulsars. We 
argue that subpulses are emitted within narrow flux tubes of the dipolar fleld 
lines and that the mean pulsar beam has a conal structure. In such a model 
the longitudes of profile components are determined by the intersection of the 
line of sight trajectory with the frequency-dependent cones of maximum mean 
intensity, while the longitudes of subpulses are determined by the intersection 
of the line of sight trajectory with subpulse-associated emission beams. Thus, 
we show that the difference in the frequency dependence of subpulse and profile 
component longitudes is a natural property of the conal model of pulsar emission 
beam. We support our conclusions by numerical modelling of pulsar emission, 
using the known parameters for this pulsar, which produce results that agree 
very well with our dual frequency observations. 

Subject headings: pulsar: subpulses: mean profiles: 

1. Introduction 

It is generally believed that the frequency dependence of pulsar radiation patterns 
(subpulses, profile components etc.) reveals the so called radius-to-frequency mapping 
effect : r(z/) oc where r is the emission altitude corresponding to the frequency u and 
p is a positive exponent. The coherent radio emission from pulsars must be generated 
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due to some instability occuring within the stream of relativistic plasma flowing along 
dipolar field lines (e.g. Kazbegi et al. 1996; Asseo & Melikidze 1998; Melikidze et al. 
2000). It is commonly assumed that the emission process is narrow-band, that is, at a 
given radio emission altitude r a relatively narrow band of frequencies, Ai/, centered at a 
certain (> Ai/ ) is generated (Cordes 1978, 1992; Kijak & Gil 1998). As a result of the 
radius-to-frequency mapping and the diverging nature of dipolar field lines, the longitudes 
of details of the observed emission (such as peaks of subpulses and profile components, 
profile edges etc.) vary with frequency. The simplest result of this is the observation that 
the overall pulse width is usually braoder at lower frequencies. One should realize that the 
actual radius-to-frequency mapping can be revealed only if the radiation processes can be 
traced along approximately the same dipolar field line (or narrow bundle of adjacent field 
lines). For this, one has to locate a narrow feature in the radiation pattern (pulse window), 
which is apparently correlated in terms of intensity at different observing frequencies. The 
frequency dependence of the longitude of occurrence of such correlated emission features is 
given by </?(i^) oc p oc r^^^{i') oc where p is the opening angle of dipolar field lines. 

The longitudinal phase ^p{iy) = 360° ■ t{v)/ P is determined by the time of arrival t{v) 
of a particular detail of emission pattern within the 360° pulse window, where P is the 
pulsar period. Usually the reference time of arrival (to = 0) corresponds to the fiducial 
plane containing the pulsar spin axis, f2, the magnetic axis, m, as well as the observer, 
O (Fig. 1). The time of arrival of radiation corresponding to a particular feature emitted 
near the fiducial plane does not depend on frequency (if the dispersive delays are properly 
removed. Craft 1970). In this paper we will use the term "detail of emission pattern" in 
several different meanings, such as: the longitude of subpulse peak ipsi the longitude of the 
mean profile peak 9?^, the longitude of the profile midpoint (prm and mean profile width 
W{i') — 2 • ipw{i^) at some intensity level w referred to the maximum intensity (usually 10% 
or 50%), where ^Pwi''^) is the longitudinal phase of mean profile corresponding to the w 
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percent of relative intensity. 

Izvekova et al. (1993) were the first to examine the frequency dependence of subpulse 
versus profile component variations for four pulsars: PSRs B0031— 07, B0320+39, B1133+16 
and B2016+28, which are well know for the prominent subpulse drift phenomenon (see 
Table 1). Izvekova et al. (1993) used a technique which did not require simultaneous 
single pulse observations at different frequencies. They measured the parameter P2, which 
is the spacing between two bands of drifting subpulses, for different frequencies and 
converted it into a longitude separation, Ay^^ = 360° ■ P2/ P- In the corresponding mean 
profiles they measured the pulse widths A(^5o = 2(^5o and A(^io = 2(^io at 50 and 10 
per cent of the maximum intensity, respectively. Izvekova et al. (1993) found that if one 
presents the frequency dependence of pulsar radiation patterns in the exponential form 
^(Pw oc then p^Q > pio > Ps- In other words, they demonstrated that the longitudes 

ips of subpulse peaks vary with frequency much less than the longitudes of corresponding 
profile components (^50 or ipiQ. Prom their sample of four pulsars, typcial values obtained 
are ~ 0.35 and Ps ~ 0.1 (Table 1). 

As Izvekova et al. (1993) noticed correctly, their results imply that the excitations 
of subpulse emission do not follow exactly the distribution of excitations leading to the 
mean emission, which is difficult to reconcile with a simple pulsar picture. However, 
Gil & Krawczyk (1996) demonstrated that this is exactly what one should expect if the 
mean pulsar beams have a conal structure. In fact, within a conal model of pulsar beams 
(Rankin 1993; Gil et al. 1993; Gil & Krawczyk 1996; Mitra & Deshpande 1999) the subpulse 
enhancement follows a narrow bundle of dipolar field lines, while the mean emission 
enhancement (corresponding to peaks of the mean profile components) is distributed over a 
cone (or a number of nested cones) of dipolar field lines (sec Fig. 1 for illustration). Thus, 
the frequency dependence of subpulse and component peak longitudes should, in general. 
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be different. On the other hand, for a patchy model of the mean pulsar beam (Lyne & 
Manchester 1988; Manchester 1995), both subpulse enhancements and those related to the 
profile components follow more or less the same bundles of dipolar field lines. Within such 
a model, the frequency dependence of subpulse and component peak longitudes should be 
similar, which is in conflict with available observational data. 

In this paper we further discuss the problem of frequency dependence of pulsar radiation 
patterns, analyzing new single pulse data obtained for PSR B0329+54 simultaneously at 
238 and 610 MHz, using the Giant Metrewave Radio Telescope (GMRT). We confirm the 
conclusion obtained by Gil & Krawczyk (1996) who used averaged data in their analysis, by 
showing the corresponding effects in single pulse data. 

2. Single pulse simultaneous dual frequency data of PSR B0329+54 

The single pulse simultaneous dual frequency data used in our analysis were taken 
with the Giant Metrewave Radio Telescope (GMRT) at Khodad in India (see Swarup et al. 
(1990) for details about the GMRT). The data were obtained by adding the total power 
detected signals from 12 of the 30 antennas of the GMRT, using the incoherent mode array 
combiner of the GMRT (see Gupta et al. (2000) for more details about the pulsar mode 
of operation of the GMRT). The bandwidth used was 16 MHz, divided into 256 spectral 
channels by the digital back-ends. The raw data were integrated to a time resolution of 
0.516 miUisecond (corresponding to a pulse longitude resolution of 0.26 degree) before 
being recorded for off-line analysis, where the data were dedispersed and gated to obtain 
the single pulse data sequence. A special feature of these observations was that they were 
carried out at two radio frequencies simultaneously, by splitting the array of 12 antennas in 
to two sub-arrays - one consisting of 5 antennas operating in the 238 MHz band and the 
other consisting of 7 antennas operating in the 610 MHz band. After total power detection. 
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the data from antennas at both frequencies were added together in the incoherent mode 
array combiner. During the off-hne analysis, the pulsar's dispersion curve was used to 
discriminate between the longitude location of the pulses at the two radio frequencies, thus 
allowing the recovery of dedispersed pulse trains at the two different radio frequencies from 
a single raw data set. The technique works as long as the combination of the two chosen 
frequencies, and the DM and period of the pulsar is such that pulses from the two radio 
bands (after dispersion) do not overlap. The advantage of the technique is that the data 
from the two frequencies is completely aligned in time, except for the dispersion delay 
between the two radio frequencies. This delay is easily and accurately estimated and the 
two dedispersed data trains are aligned by removing this delay. The data set used here was 
recorded on 10th August, 1999 and consists of a sequence of 1475 pulses of PSR B0329+54. 

The average profiles for PSR B0329+54 obtained at 238 and 610 MHz from this data 
are shown in Fig. 2. Both profiles show a central core component and the two prominent 
conal components (one on each side of the core) that are normally easily seen for this 
pulsar. The radius-to-frequency mapping effect for the average profiles is also clearly seen 
in this figure, with the conal components at the lower frequency (238 MHz) occuring further 
away from the central core component than the corresponding conal components at the 
higher frequency (610 MHz). A sample of selected single pulses at the two frequencies is 
shown in Fig. 3, along with the integrated profiles for reference. The data for 238 MHz is 
shown by the thicker hues, whereas the 610 MHz data is shown by the thinner hues. These 
single pulses are selected on the basis of good signal to noise ratio (at both frequencies) 
for at least one of the conal components. The threshold used was a peak SNR criteria (in 
the longitude window for the conal component) of 7 times the rms of the off-pulse noise. 
This yielded 20 and 16 pulses for the leading and trailing components respectively, from 
which the 5 best examples - where the sub-pulse shapes at both frequencies are fairly clear 
and Gaussian-like - were selected for each conal component. The longitude windows used 
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for the two conal components are shown in the figure by the grey vertical fines in tfie left 
fiand panels, and are demarcated as "winl" and "win2" for tfie leading and trailing conal 
components, respectively. 

It is seen tfiat tfie emission features in tfie individual pulses also exfiibit tfie radius- 
to-frequency mapping effect in tfiat tfie subpulses at tfie lower frequency appear earlier for 
tfie leading conal component and later for tfie trailing component, wfien compared witfi 
tfie corresponding feature at 610 MHz. However, tfie interesting point is tfiat, for all tfie 
single pulses sfiown, tfie longitudinal separation between tfie subpulse peaks at tfie two 
frequencies is smaller tfian tfie corresponding separation in tfie average profile. Tfiis is 
quantified by examining tfie peaks of tfie cross-correlation functions (CCFs) for tfie average 
profile as well as tfie individual pulses. Tfiese CCFs are displayed in tfie rigfit fiand panel 
of Fig. 3. Tfie peaks of tfie CCFs (marked by tfie open circle symbols) are a measure of tfie 
longitudinal separation between tfie correlating single pulse features at tfie two frequencies, 
for tfie selected conal component. For tfie average profile, tfie location of tfie CCF peaks are 
marked by tfie two grey vertical lines in tfie rigfit fiand panel. Tfiese correspond to a mean 
sfiift of 7 bins (1.82 degrees) of pulse pfiase (negative) for tfie leading conal component 
and 4 bins (1.04 degree) of pulse pfiase (positive) for tfie trailing conal component, for 
tfie average profile. Tfiese two sfiifts are not equal in amplitude, witfi tfie departure from 
symmetry being about (7 — 4)/2 = 1.5 bins or 0.39 degree. Tfiis asymmetry, wfiicfi is 
also clearly visible to tfie eye in Fig. 2, is believed to be due to tfie effects of retardation 
and aberration (see Gangadfiara & Gupta 2001) wfiicfi cause tfie conal beams at different 
frequencies (emitted at different fieigfits) to be aberrated and retarded by different amounts, 
leading to increased separation for tfie leading parts of tfie profile and reduced separation 
for tfie traifing parts of tfie profile. Tfie retardation-aberration sfiift of 0.39 degree or 0.77 
msec, translates into a difference of emission altitudes of Ar ~ c7.7 x 10~^s ~ 2.3 x 10^ cm, 
wfiere c = 3 x 10^° cm/s. Tfiis is perfectly consistent witfi tfie empirical radius-to-frequency 
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mapping given by Kijak & Gil (1998). In fact, for PSR B0329+54 Kijak & Gil (1998) give 
the emission altitude r(z/) = (5 ± 0.5) x i^ghI^'^'^'^'^ cm. For our frequencies of 0.24 and 
0.61 GHz, the emission altitudes are about 7 x 10^ cm and 5 x 10^ cm, respectively. Thus, 
the difference in altitudes is Ar ~ 2 x 10^ cm, in good agreement with Ar ~ 2.3 x 10^ cm 
inferred from the observed retardation-abberation shift. These numbers are also in good 
agreement with the estimates for emission heights for the conal components for this pulsar 
obtained by Gangadhara & Gupta (2001) using accurate measurements of the retardation 
and aberration effects. 

The typical deviation of the peak of the CCF for the individual pulses with respect 
to that for the average profile for the corresponding conal component, is 2 to 3 bins (0.52 
to 0.78 degree) of longitude, with a minimum to maximum range of to 5 bins (0 to 
1.3 degree). There appears to be some evidence that these shifts in the CCF peaks are 
systematically larger for the trailing conal component. Note that there are practically no 
cases where the magnitude of the shift for the individual subpulse is more than that for the 
corresponding conal component in the average profile. 

The above conclusions are further quantified by studying the distribution of the CCF 
peaks for the individual pulse cases, compared to the peak of the CCF for the average 
profile for the corresponding conal component. This is illustrated in Fig. 4 which shows the 
histograms of the distributions of the CCF peaks obtained for the selected single pulses for a 
SNR threshold of 5 times the off-pulse rms noise. The results for the two conal components 
are shown as separate histograms, for the respective windows. The total number of pulses 
that meet the threshold criteria are 56 and 57 for the two conal components, out of a total 
of 1475 pulses. The actual quantity used for the histogram plots is the difference between 
the CCF peak for the individual pulse and the CCF peak for the average profile for the 
corresponding conal component. Thus the zero position bin in the histogram corresponds 
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to cases of single pulses where the shift in the subpulse at the two frequencies is the same 
as the corresponding shift in the average profile. These plots clearly show that the CCF 
peaks for the individual pulses are asymmetrically distributed with respect to the CCF 
peak for the average profile for the corresponding conal component. They also show that 
the magnitude of the shift is somewhat larger for the trailing component, as compared to 
the leading component, with mean values for the shifts being -2.5 bins (-0.65 deg) and 1.4 
bins (0.36 degree), respectively. 

These results are consistent with those of Izvekova et al. (1993). All these are 
summarized in Table 1. Although there are differences from pulsar to pulsar, the results 
in Table 1 reinforce the conclusion that the frequency dependence of subpulse patterns is 
apparently weaker than the frequency dependence of corresponding features in the average 
profile. The exponents quoted for PSR B0329-I-54 (0.14 and 0.07 for mean profile and single 
pulses, respectively) refer to average values following from estimates obtained for the left 
(L) and right (R) side of the pulse profiles, after corrections for retardation-aberation shift 
of about 1.5 bin (0.4 degree of longitude). In the next section we argue that this different 
frequency dependence of subpulses and corresponding profile components is natural, 
provided that the mean pulsar beam has a conal structure and the subpulse emission 
is associated with narrow bundles of dipolar magnetic field lines. We assume that only 
one bundle (controlling the subpulse-associated plasma flow) is related with one subpulse 
observed simultaneously at the two frequencies. 



3. Structure of pulsar beams 

It is widely believed that pulsar elementary emission is relativistically beamed tangently 
to dipolar magnetic field lines. Thus, the geometry of the emission region can be described 
by the opening angle p(i/) = 1°.25 s r6(i^)^^^P~^^^ (Fig- 1); where P is the pulsar period. 
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and re = r{i')/R is the frequency dependent emission altitude in units of stellar radius 
R = 10^ cm. The mapping parameter (0 < s < 1) is determined by the locus of dipolar 
field lines on the polar cap (s = at the pole and s = 1 at the polar cap edge). Thus, 
s — d/vp, where d is the distance from the magnetic axis, m, to a field line on the polar cap 
corresponding to a certain detail of the emission pattern, and Vp — 1.4 • 10^P~^/^ cm is the 
canonical polar cap radius. According to the widely accepted concept of radius-to-frequency 
mapping, the higher frequencies are emitted at lower altitudes than the lower frequencies. 
Kijak & Gil (1998) have found an empirical formula describing an altitude of the emission 
region corresponding to a frequency z/ghz (in GHz) as 

re ^ 50 ■ .^°f . re-°-°^ • P' '' (1) 

where Tq is the pulsar characteristic age in million years and P is the pulsar period. This 
empirical relationship will be used in our model calculations. 

To perform geometrical calculations of the radiation pattern one has to adopt a model 
of instantaneous energy distribution on the polar cap. We will assume that at any instant 
the polar cap is populated by a number of features with a characteristic dimension V (e.g. 
Gil & Sendyk 2000) , each delivering to the magnetosphere a column of plasma flowing along 
separate bundle of dipolar magnetic field lines. These plasma columns are supposed to be 
sources of the narrow-band radiation, with different radio frequencies emitted at different 
altitudes, as described by the radius-to-frequency mapping law (Kijak & Gil 1998). Thus, 
one plasma column illuminates different spots at different frequencies. Some of these spots 
can be intersected by the line of sight (LOS henceforth) and thus observed simultaneously 
at two frequencies. Each feature is modelled at the polar cap by a gaussian intensity 
distribution 

h = exp(-«/Vp2) (2) 
where — d'^{ip) + df — 2d{ip)di ■ cos[a{(p) — ai]. Here the pairs [dj, ai] and [d (</?), cr (</?)] are 
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polar coordinates of the i-th feature (maximum intensity) and the observer's position on the 
hne of sight (projected onto the polar cap), respectively. The instantaneous emission of the 
n-th pulse is described by /„((^) = E^'\/j((^), where is the number of adjacent features 
contributing to the radiation observed in the n-th pulse at a longitude (f. The average pulse 
profile is therefore I{ip) — ■^T,^=iln{<f): where A^" is the number of averaged single pulses. 
The azimuthal angle variations of the i-th feature are described by 

(jj = arctan[sin ipi ■ sm{a + /?) ■ sin a/ cos(q; + /?) — cos a ■ cos pi] — n ■ Act , (3) 

where 

p = 2 ■ asin[sin2((^/2) ■ sina • sin(Q; + (3) + sin2(/?/2)]^/^ = l.°25 s rQ{vf''^p-'^''^ , (4) 

and Act is the change of azimuthal angle per one pulsar period, a is the inclination angle 

between the rotation and magnetic axes and (3 is the impact angle - the angle of the closest 
approach of the line of sight to the magnetic axis (Figs. 1, 5 and 6). 

Any specific intensity distribution can be transferred along dipolar field lines 
to the emission region and then along straight lines to a given observer specified 
by inclination and impact angles {a, (5). Thus, knowing a distribution of opening 
angles Pi{i') — l.°25((ij/rp)rg^^(z/)P~-'^/^ corresponding to instantenous radiation, 
one can calculate the longitude of an i-th detail of the emission pattern ipi = 
2 asin { [sin (^^) • sin (^^)] /[sina ■ sin(Q; + /9)]}^^^ at a given frequency v (Gil 1981). 
The mean profile can be obtained by a simple averaging of individual waveforms, which can 
in principle vary in time, depending on the adopted model of distribution of instantaneous 
arrangement of energy on the polar cap. 

Gil & Krawczyk (1996) and Gil & Sendyk (2000) argued that the nested cone structure 
of mean pulsar beams is induced by the vacuum gap E x B drift (Ruderman & Sutherland 
1975), which forces the subpulse-associated plasma columns to rotate around the magnetic 
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axis. This model suggests a simple and natural relationship between a sequence of single 
pulses and features of mean emission. It is also very convenient for geometrical simulations, 
because the azimuthal angle cr (eq. 3) of each plasma column varies regularly with time 
(sequential pulse number). The frequency dependence of pulsar emission patterns within 
such a model is qualitatively illustrated in Fig. 1. The subpulse-associated plasma columns 
corresponding to the half power beam width of the subpulse emission (thin small circles 
marked by Slu)), which will be called spots, perform a circumferential motion around 
the magnetic axis m, in which the azimuthal angle a varies with time while the opening 
angle p(z/) remains constant. This motion determines the cones (thick large circles) of 
the maximum average intensity with the opening angles pi and p2 at frequencies ui and 
1/2, respectively. The frequency dependent longitude (fpiv) of the profile component is 
determined by the intersection of the line of sight with the average cone at the frequency 
dependent opening angle p(i^). The frequency dependent longitude v^s(z^) of subpulse peak 
is determined by the local maximum intensity along the cut of the line of sight through 
the subpulse spot. These cuts are marked by dashed lines (shorter dash for S'{i'i) and 
longer dash for S{i'2), respectively). The position of the spot S was choosen in such a way 
that (fisi^i) = ^pi^i)- This is however a special case, which does not hold for the spot 
5", in which case the longitude of the subpulse peak (p's{i^i) is different from the longitude 
of a corresponding component ipp{i'i). When the frequency changes from ui to ^2, the 
subpulse peak longitude changes from ips{i^i) to (/9s(z/2), while the corresponding profile peak 
longitude changes from (pp{i'i) to (pp{u2)- Note that A(^p — (pp{v2) — ^p{yi) is generally 
larger than Ac/?^ = Vs{''^2) ~ </'s(^i)- This is consistent with observations of four pulsars 
pubhshed by Izvekova et al. (1993) and confirmed in this paper for yet another pulsar PSR 
B0329+54 (Table 1). 

It is important to note that in this model, ^(^s < Ay^p, is always true. The equality 
condition holds only for two special cases: (i) a ~ 90° and /3 ~ 0° and (ii) cr ~ 0°. The 
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latter case is trivial since it corresponds to Ac/?^ = A(/7p ^ 0° for which both subpulses are 
observed at the fiducial phase; while the former one corresponds to central cut of the LOS 
across the beam. For all cases other than these two, Aips < A(^p is strictly true. This 
prediction of the model has been shown to be true in the observed data (Figs. 2 and 3) and 
will be shown to hold for the simulated synthetic data as well. It should be emphasized that 
we consider only cases in which both higher and lower frequency spots illuminated by the 
same single subpulse-associated plasma column are observed (i.e. intersected by the LOS). 
The cases in which subpulse at one of the two frequencies is missing are not interesting for 
the present analysis, as no frequency dependent effect can be measured for these. 

After qualitative illustration of the problem (Fig. 1) we now explore it quantitatively 
by means of geometrical simulations. Although we did not intend to model exactly the 
pulse profiles of PSR B0329+54, we adopted in our calculation the values of P = 0.7 s 
and P — 2 ■ 10~^^. Thus, our model of the polar cap or the corresponding instantaneous 
radiation beam presented in Fig. 5 corresponds in a sense to PSR B0329+54. We followed 
Gil & Sendyk (2000) (their eqs. [11] and [12]) for description of the polar cap arrangement. 
The characteristics of instantaneous radiation are imprinted on the polar cap, transmitted 
along dipolar field lines to the radio emission altitudes and then emitted to a particular 
observer along straight paths. The circumferential motion of sparks around the magnetic 
axis (Fig. 5) will lead to a nested cone structure of the average radiation beam, as presented 
in Fig. 6. The adopted inclination angle a — 30° and four line of sight trajectories 
correspond to impact angles P — 0, 3.3, 5.0 and 6.6 degrees, respectively. We performed 
calculations for the two frequencies 238 and 610 MHz. Hence, each line of sight trajectory 
bifurcates when projected onto the surface of the polar cap. 

The results of calculations of radiation waveforms corresponding to emission patterns 
presented in Figs. 5 and 6 are shown in Fig. 7. Four vertical panels correspond to different 
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impact angles (3 — 0, 3.3, 5.0 and 6.6 degrees, respectively. The broader emission at 238 
MHz is shown by dashed lines and higher frequency radiation at 610 MHz is shown by solid 
lines. Both instantaneous single pulse and mean profile are presented for each impact angle. 
We have calculated a long sequence of single pulses to make sure the mean profile is stable. 
We selected pulse #4 in our sequence (for no particular reason) and show its profile for 
each impact angle at the two frequencies. 

The first panel corresponds to a central cut of the line of sight with the impact angle 
= 0° and, of course, a ~ 90°. In this case the longitudinal difference of subpulse and profile 
component peaks occuring at different frequencies is exactly the same [Aip^^^ — Aipp). 
This can be seen from both visual inspection of the pulse profiles and from the CCF 
analysis. The CCF computed separately for left (L) and right (R) component, peaks at 
the same lag in both mean and single pulses. However, as (3 increases, the difference for 
subpulses becomes smaller and smaller as compared with a corresponding difference for 
peaks of profile components. This translates into CCF lags which are smaller for single 
pulses than for mean pulses. The lag difference is about 0.5, 1.0 and 2.0 degrees for (3 values 
of 3.3, 5.0 and 6.6 degrees, respectively. This can be compared with Fig. 3, where the mean 
lag difference is about 2 sample bins or 0.5 degrees of longitude. Thus, the actual data 
presented in Fig. 3 correspond to a modelling case with /3 ~ 3 degrees (see Fig. 7). 

Another interesting, though obvious property is that the frequency separation 
disappears near the fiducial longitude = 0°, as in the actual data presented in Figs. 2 
and 3. We have selected three central subpulses marked by 1, 2 and 3 in Fig. 7. The CCF 
computed for these subpulses show lags near zero in Fig. 8 (for the subpulse 1 the lag is 
exactly zero) . The explanation of this effect is simple and natural. The fiducial logitude was 
chosen in such a way that it corresponds to the fiducial plane containing both the ft and 
m axes as well as the observer O. The divergence of dipolar magnetic field lines does not 
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affect longitudinal phases of the observed radiation emitted exactly at this plane (subpulse 
1) and affects it only slightly near this plane (subpulses 2 and 3). 

In our simulations we have not included the effects of retardation and aberation, 
which were unequivocally detected in our high resolution, simultaneous dual frequency 
single pulse data. These effects shift the entire low frequency emission patterns towards 
earlier phases (leading side) with respect to high frequency patterns. This results in larger 
frequency difference on the left (L) than on the right (R) side of the profile, thus producing 
the asymmetrical CCF lags seen in Fig. 3. However, the deviation of single pulse lags from 
mean pulse lags does not depend on the retardation-aberation shift (see Fig. 2), so our 
symmetrical model CCF can be compared with the asymmetrical CCF computed for the 
actual data. 



4. Discussion and conclusions 

The observed emission patterns of pulsar radio emission are usually frequency 
dependent. The longitudes of details of pulse profiles appear earlier/later for the 
leading/trailing side of the pulse window, as frequency decreases. This is well understood 
in terms of the radius-to-frequency mapping and diverging nature of dipolar magnetic field 
lines in the radio emission region. However, in PSR B0329-I-54 and four other pulsars 
(Table 1) it has been observed that for emission at two different frequencies, the difference 
in longitudes of subpulse peaks is generally smaller than the difference in longitudes of 
corresponding component peaks. This suggests that the subpulse and profile component 
excitations do not follow the same bundles of dipolar field lines, which appears surprising at 
first sight. However, we argued in this paper that this is exactly what one should expect if: 

(i) the subpulse enhancements are related to relativistic plasma columns following a 
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narrow (and separated from each other) bundle of dipolar field lines, with different 
frequencies generated at different altitudes rQ{i') = r{u)/R, 

(ii) the profile component enhancements are related to the conal structure of the mean 

pulsar beam, and 

(iii) the hne of sight trajectory of the observer across the beam is neither central not 
grazing (i.e. the impact angles < /3 < p;,, where = 1°.25 rg'^^P"^/^ is the opening 
angle of the full pulsar beam). Particularly, for PSR B0329+54 we deduced the 
impact angle /3 ~ 3°, which agrees well with other independent estimates from about 
2° (Rankin 1993) to about 4° (Gil and Lyne 1995). 

In deriving our conclusions, we have made a number of well justified assumptions, 
which we summarize below: 

(1) The pulsar radiation is relativistically beamed along dipolar lines of a global magnetic 

field of the neutron star. This is a generally well accepted assumption. 

(2) The narrow-band pulsar radio emission follows the radius-to-frequency mapping 

r — r{y). In our calculations we used a specific form of the empirical relation obtained 
by Kijak & Gil (1998). Although, the actual form of the radius-to-frequency mapping 
is not important for our general conclusions (as long as the higher frequencies are 
generated closer to the polar cap than the lower ones), this relation is strongly 
supported by the retardation-aberation shift clearly revealed in our high-resolution, 
simultaneous dual frequency data. 

(3) The instantaneous emission beam consist of a number of subpulse-associated beams 

emitted by the spark-associated plasma columns. Thus, we invoke a non-stationary 
model of inner accelerator discharge via a number of isolated sparks (Sturrock 1971; 
Ruderman & Sutherland 1975; Gil & Mitra 2001). 
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(4) The structure of mean pulsar beam results from a behaviour of subpulse- associated 
beams averaged over time intervals much longer than the pulsar period. In our model 
these beams perform a circumferential motion around the magnetic axis following the 
E X B drift motion of sparks on the polar cap. Strong observational evidence of such 
circumferential motions has been recently presented by Deshpande & Rankin (1999, 
2001). 

The model of the polar cap arrangement we have adopted here is highly symmetric, 
which simplifies geometrical calculations. However, our general conclusions do not depend 
on any specific details of this model, as long as the instantaneous beam is organized into a 
group of isolated subbeams (Fig. 5) and the average beam has a conal structure (Fig. 6). In 
our model this conal structure is caused by the circumferential motion of spark-associated, 
subpulse- producing plasma columns around the axis of symmetry (magnetic axis). For 
a convenience, we have calculated the emission patterns using a specific velocity of 
circumferential motion following from the E x B drift applied to sparks (eq. [14] in Gil & 
Sendyk 2000), although the value of this velocity is not critical for our results. Our general 
conclusions apply to a more general picture of the organization of pulsar beams, in which 
enhancement of excitations leading to subpulses and corresponding profile component, 
respectively, do not follow each other. We consider the E x B drift induced model as the 
simplest and most natural, but the reason for the formation of the conal structure of the 
mean pulsar beam is not important for our conclusions. We would also like to emphasize 
that our conclusions do depend sensitively on the impact angles (5 of the closest approach 
of the fine of sight to the magnetic axis m (Fig. 1) but they are rather independent of the 
inclination angles a. 

To conclude, we have demonstrated in this paper that the angular beaming model 
of instantaneous emission accompanied by the conal model of mean emission predict 
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frequency dependent properties of pulsar radiation that match well with those observed in 
the simultaneous, dual frequency single pulse data of PSR B0329+54 that is reported here. 
This model also explains the earlier observations reported for four other pulsars - PSRs 
B0031-07, B0320+39, B1133+16 and B2016+28. We therefore conclude that, at least 
in these pulsars, the mean pulsar beams have a conal structure, possibly induced by the 
circumferential E x B drift motion of subpulse-associated beams. Thus, our paper provides 
independent support for the nested-cone model of pulsar beams (Rankin 1993; Gil et al. 
1993; Kramer et al. 1994; Mitra & Deshpande 1999; Gangadhara & Gupta 2001). 

This paper is supported in part by the Grant 2 P03D 008 19 of the Polish State 
Committee for Scientific Research. We thank an anonymous referee for useful comments. 
We are thankful to Drs. D. Mitra and B. Zhang for helpful discussions. We thank E. Gil, 
O. Maron and M. Sendyk for technical assistance. We also thank the staff of the GMRT 
that made these observations possible and specifically Dr. B.C. Joshi for this particular 
data set. GMRT is run by the National Centre for Radio Astrophysics of the Tata Institute 
of Fundamental Research. 



- 19- 

REFERENCES 

Asseo, E., Melikidze, G. 1998, MNRAS, 301, 59 

Cordes, J.M. 1992, in JAY Coll. 128, ed. T.H. Hankins, J.M. Rankin & J.A. Gil (Zielona 
Gora: Pedagogical Univ. Press), 253 

Cordes, J.M. 1978, ApJ, 222, 1006 

Craft, H.D. 1970, Kept. No. 399, Cornell University 

Deshpande, A.A., & Rankin, J.M. 1999, ApJ, 524, 1008 

Deshpande, A.A., & Rankin, J.M. 2001, MNRAS, 322, 438 

Gangadhara, R.T.G., & Gupta Y. 2001, ApJ, 555, 31 

Gil, J. 1981, Acta Phys. Pol, B12, 1081 

Gil, J., Kijak, J., Seiradakis, J.H. 1993, A&A, 272, 268 

Gil, J., Lyne, A.G. 1995, MNRAS, 276, L55 

Gil, J., & Krawczyk, A. 1996, MNRAS, 280, 143 

Gil, J., Sendyk, M. 2000, ApJ, 541, 351 

Gil, J., & Mitra, D. 2001, ApJ, 550, 383 

Gupta, Y., Gothoskar, P. G., Joshi, B. C, Vivekanand, M., Swain, R., Sirothia, S., & Bhat, 
N.D.R. 2000, in lAU Colloquium 177: Pulsar Astronomy - 2000 and Beyond, eds. 
M. Kramer, N. Wex & R. Wielebinski, ASP Conf. Series, 105, 277 

Izvekova, V.A., Kuzmin, A.D., Lyne, A.G., Shitov, Yu.R, & Smith, F.G. 1993, MNRAS, 
261, 865 



-20- 

Kazbegi, A., Machabeli, G., Melikidze, G., Shukre, C. 1996, A&A, 309, 515 
Kijak, J., Gil, J. 1998, MNRAS, 299, 855 

Kramer, M., Wielebinski, R., Jessner, A., Gil, J. A., Seiradakis, J.H. 1994, A&ASuppl. Ser., 
107, 515 

Lyne, A.G., Manchester, R.N. 1998, MNRAS, 234, 477 
Manchester, R.N. 1995, JA&A, 23, 283 

Melikidze, G.I, Gil, J., & Pataraya, A.D. 2000, ApJ, 544, 1081 

Mitra, D., & Deshpande, A.A. 1999, A&A, 346, 906 

Rankin, J.M. 1983, ApJ, 274, 333 

Rankin, J.M. 1990, ApJ, 352, 247 

Rankin, J.M. 1993, ApJ, 405, 285 

Ruderman, M.A., Sutherland, RG. 1975, ApJ, 196, 51 

Sturrock, RA. 1971, ApJ, 164, 529 

Swarup, G., Ananthakrishnan, S., Subrahmanya, C. R., Rao, A. P., Kulkarni, V. K., & 
Kapahi, V. K. 1990, in High Sensitivity Radio Astronomy, ed. N. Jackson & R.J. 
Davis (Cambridge Univ. Press) 

Usov, V.V. 1987, ApJ, 320, 333 



This manuscript was prepared with the AAS macros v5.0. 



- 21 - 



Table 1. Frequency dependence of pulse profile characteristics ipi (xv ^\ 



PSR B 


(^50 OC P 








Ref. 


0329+54 






0.14 


0.07 


(2) 


0031-07 


0.45 


0.35 




0.05 


(1) 


0320+39 


0.34 


0.18 




0.15 


(1) 


1133+16 






0.25 


0.06 


(1) 


2016+28 


0.35 


0.2 




0.17 


(1) 



Note. — (fi5Q - half width W^q — 2 ■ (f^o of the profile at 50% maximum 
intensity, ipio - half width Wiq = 2 ■ cpiQ of the profile at 10% maximum 
intensity, ipc - longitude of profile component peak, ipg - longitude of 
subpulse peak, Ref: (1) Izvekova et al. (1993), (2) this paper 
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Figure captions 

Fig. 1. — The conal model for the mean pulsar beam. The geometry of observation is 
determined by the inclination angle a between the magnetic axis m and the spin axis Q, and 
the impact angle (3 of the closest approach of the line of sight (observer O) to the magnetic 
axis m. The plane containing the spin axis O, the magnetic axis m and the observer O 
defines the fiducial longitude ip — Q° . The frequency dependent position of any observed 
feature in the beam pattern is described by two angles: the frequency dependent opening 
angle piy) between the m axis and the line of sight and by the azimuthal angle a between 
the fiducial plane and the plane of dipolar field lines associated with a particular feature. 
The subpulse emission is associated with the frequency dependent subpulse spots S{y) and 
the profile components are associated with the frequency dependent cones. The longitudes 
of subpulse peaks (ps{^) and profile peaks ipp^u) are marked, with U2 > ui. 

Fig. 2.— Average profiles at 238 and 610 MHz for PSR B0329+54, obtained from the 
simultaneous dual frequency observations. The darker curve is for 238 MHz and the lighter 
one for 610 MHz. Each sample bin corresponds to 0.26 degrees of longitude. 

Fig. 3. — The observed frequency dependence of radio emission patterns in PSR B0329+54. 
Selection of 10 single pulses from the observations is shown (at both frequencies) in the 
left hand panels. The darker curves are for 238 MHz and the lighter curves for 610 MHz. 
The average profile (top panel) is shown truncated in amplitude to highlight the conal 
components. The right hand panels show the cross-correlation functions for the relevant 
conal component, with the peaks marked by the open circle symbols. Each bin or lag 
corresponds to 0.26 degrees of longitude. 

Fig. 4. — Histograms of the distributions of the location of the CCF peaks obtained for the 
leading conal component (top panel) and trailing conal component (bottom panel) for the 
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selected single pulses. The locations are taken with respect to the location of the peak of 
the CCF for the average profile for the corresponding conal component. 

Fig. 5. — Instantaneous arrangement of sparks on the Polar Cap (or instantaneous pulsar 
beam projected onto the PC along dipolar field lines) corresponding to the simulated pulse 
#4 presented in Fig. 7. Four different (frequency bifurcated) line of sight trajectories with 
impact angles (3 — 0.0, 3.3, 5.0, 6.6 degrees corresponding to frequencies 238 and 630 MHz 
are presented. The pulsar rotates clock-wise, so the leading part of the pulse corresponds to 
the right side of this figure (or left if viewed from the spin axis il) . The inclination angle 
between the magnetic axis (marked by the central dot) and the spin axis is 30 degrees. 

Fig. 6. — The average structure of the polar cap (or mean pulsar beam projected onto 
the PC) resulting from the circumferential motion of sparks (or subpulse beams) around 
the magnetic axis (marked by the central dot). For explanation of different line of sight 
trajectories see caption for Fig. 5. 

Fig. 7. — The frequency dependence of pulsar radiation patterns corresponding to Figs. 5 
and 6. In the analogy to Fig. 3, the left panels present the pulses, while the right panels 
present the corresponding cross-correlation functions (CCF). Each vertical panel presents 
the integrated profile and a single pulse for a given impact angle 13 = 0.0, 3.3, 5.0 and 6.6 
degrees. The symbols L and R correspond to the leading and trailing conal components 
of the profile, respectively. One can easily notice from both visual inspection and from 
the CCF analysis that the difference in longitudes of pulse features occuring at 238 and 610 
MHz depends on the impact angle f3 (closest approach of the observer to the beam axis). For 
(3 = 0, this separation is the same for both single pulses and integrated profiles. However, at 
larger impact angles the separation for subpulses is smaller than for a corresponding profile 
components. It is worth noting that longitudes of subpulse peaks and corresponding profile 
components are the same near the fiducial phase ip = 0° like in the real data presented in 
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Figs. 2 and 3 (see CCF analysis in Fig. 8 corresponding to central subpulses marked by 1, 
2 and 3 in this figure). 

Fig. 8. — The cross-correlation function (CCF) calculated for three subpulses marked by 1, 
2 and 3 in Fig. 7. The maximum of each function is at zero lag, meaning that subpulses 
near the fiducial longitude 0° (profile centre) appear in the same phases at both frequencies, 
unhke the side subpulses analysed in Fig. 7. 



